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Abstract. We review a coarse-graining strategy (multiblob approach) for polymer 
solutions in which groups of monomers are mapped onto a single atom (a blob) and 
effective blob-blob interactions are obtained by requiring the coarse-grained model to 
reproduce some coarse-grained features of the zero-density isolated-chain structure. 
By tuning the level of coarse graining, i.e. the number of monomers to be mapped 
onto a single blob, the model should be adequate to explore the semidilute regime 
above the collapse transition, since in this case the monomer density is very small 
if chains are long enough. The implementation of these ideas has been previously 
based on a transferability hypothesis, which was not completely tested against full- 
monomer results (Pierleoni et ai, J. Chem. Phys., 2007, 127, 171102). We study 
different models proposed in the past and we compare their predictions to full-monomer 
results for the chain structure and the thermodynamics in the range of polymer volume 
fractions $ between and 8. We find that the transferability assumption has a 
limited predictive power if a thermodynamically consistent model is required. We 
introduce a new tetramer model parametrized in such a way to reproduce not only 
zero-density intramolecular and intermolecular two-body probabilities, but also some 
intramolecular three-body and four-body distributions. We find that such a model 
correctly predicts three-chain effects, the structure and the thermodynamics up to 
<f> ~ 2, a range considerably larger than that obtained with previous simpler models 
using zero-density potentials. Our results show the correctness of the ideas behind the 
multiblob approach but also that more work is needed to understand how to develop 
models with more effective monomers which would allow us to explore the semidilute 
regime at larger chain volume fractions. 
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1. Introduction 

Polymer solutions show a wide variety of behaviors, depending on chain length, density, 
and temperature. [U [21 [31 0] In the dilute regime the isolated chain radius of gyration R g 
is the relevant length scale and the properties of the solution can be described in terms 
of single-chain properties and of the solvent- quality parameter. The radius of gyration 
scales as R g = bL u , where b is the monomer characteristic length (Kuhn segment), 
which depends on chemical details and temperature, L is the number of monomers per 
chain, and v is a universal exponent. In the good-solvent regime[5] v = 0.587597(7), 
while v = 1/2 (with logarithmic corrections [31 E]) for ^-solvents. The semidilute regime 
is entered when chains start to overlap, but still the monomer density is small. If 
c = N/V is the polymer concentration — N is the number of chains and V the volume 
of the system under consideration — and c m = cL is the monomer concentration, the 
semidilute regime is characterized by c > c* (or equivalently by $ > 1, where $ = c/c* 
is the polymer volume fraction) and c m « 1, where c* = 3/(4ttR^) is the overlap 
concentration. Note that c m = (3/47r6 3 )$L 1 ~ 3i/ . Hence, when increasing $, increasingly 
longer polymer chains are needed to ensure the semidilute condition c m < 1. If c m is 
not small, one enters the concentrated or melt regime. 

From this discussion it appears that very long polymers are necessary to obtain a 
genuine semidilute regime over several orders of magnitude in chain density. For this 
reason, simulations of semidilute solutions of linear chains, even at the level of generic 
lattice or bead-spring models with implicit solvent, are quite expensive and have been 
limited to not too long chains and not too high densities. [6], [TJ, [HI [HI UHl EH] Moreover, 
in many complex situations, polymers only constitute one species in the solution, 
making full- monomer simulations even more difficult. In these cases a modelling at 
length scales of the order of the polymer size is often sufficient to provide the relevant 
thermodynamic and structural informations on the solution. For instance, to determine 
the phase behavior of polymer-colloid mixtures in the colloid regime, [121 HSl US] a 
detailed microscopic model of the polymers is not necessary. It is enough to use coarse- 
grained models which retain the essential thermodynamic (long-wavelength) behavior 
of the polymer solution. [16] Other examples are block copolymers for which the self- 
assembling of the chains in supramolecular aggregates of various shapes and sizes is 
ubiquitous. The description of the physical behavior of the self-assembled phases only 
requires a modelling at the mesoscopic level rather than at the microscopic (monomer) 
level. [T7J HHJ [191 120] This is the realm of self-consistent field-theoretical methods which 
have proved to be very effective to describe the physics of concentrated solutions and 
melts of homopolymers and block-copolymer blends. [211 E] 

Coarse-grained models for soft condensed matter systems have received much 
attention in the last two decades. [22] In the simplest approach one maps polymer chains 
onto point particles interacting by means of the pairwise potential of mean force between 
the centers of mass of two isolated polymers. [23l [22] This potential is of the order of 
2k B T at full overlap, [2H [25] has a limited range of the order of 3R g and is very well 
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represented by a linear combination of few Gaussian functions. [HI [EE] Such a model 
is however limited to the dilute regime $ < 1, in which many-body interactions [26J 
can be neglected. This limitation was overcome ten years ago in a seminal work, j27J IB] 
which eliminated the complexity of the many-body interactions by introducing a density- 
dependent pair potential, which is unique according to Henderson theorem [28] and 
reduces to the potential of mean force in the limit of zero density. 

This work has paved the way to the use of soft effective particles to represent 
polymer coils in complex situations such as in modelling colloid-polymer mixtures. [16] 
However, density-dependent potentials are difficult to handle. Care is needed to derive 
the correct thermodynamics [29] [30], [31] and to compute free energies and phase 
diagrams. [22] Also their use in non-homogeneous situations is cumbersome since the 
interaction should depend on the local density which is not known beforehands and 
some kind of self-consistent procedure should be developed. Furthermore, representing 
polymers as soft spherically symmetric particles is not always appropriate. For instance, 
in studying polymers adsorbed on surfaces, like polymer brushes or polymer-coated 
colloids, it is clear that the anchorage to the surface breaks the rotational symmetry 
of the chains, an effect that must be taken into account in any accurate coarse-grained 
model. A further example is in modelling solutions of A-B block copolymers which 
cannot be represented as soft particles interacting by a spherically symmetric pair 
potential. [UJ EEl [191 ED] 

In principle those limitations can be overcome by switching to a model at a lower 
level of coarse graining, i.e. by mapping a long linear polymer to a short linear chain 
of soft effective monomers (called "blobs" in the following). Such model retains some 
internal degrees of freedom, which allow more flexibility in chain geometry as necessary, 
for instance, in anisotropic systems. Moreover, in the semidilute regime, this model is 
expected to allow the use of density-independent blob-blob interactions, since the local 
density of the blobs can always be kept small by increasing the number of blobs per 
chain. Indeed, if chains of L monomers are partitioned in n effective blobs of m = L/n 
monomers each, the local concentration of the blobs is q, = cn. The blob overlap 
concentration is given by c£ = 3/(47rr 3 ), where r g is the radius of gyration of the blob. 
If we assume that r g = bm u , where b is the Kuhn length that appears in the scaling of 
the radius of gyration (this relation, though not exact, is a very good approximation), 
we obtain 

°b 3 UC ^ 1-3,, 1-0.763 



c* b 4tt6 3 m 3 " 



^n 1 - 3 " « $n-°- 763 . (1) 



Hence, for any polymer volume fraction $, since v > 1/3 above the collapsed phase, one 
can choose n so that Q,/c£ < 1, i.e., so that blobs do not overlap. In this regime, the 
size of each blob is approximately density independent, hence each blob can be replaced 
by an effective single atom of fixed size in a coarse-grained representation. Moreover, 
one expects the many-body interactions among blobs of different chains to be negligible. 
Hence, the parametrization of the intermolecular interactions among the chains in terms 
of zero-density pair potentials should be reasonably accurate. Conversely, one can use 
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the n blob model with zero-density intermolecular pair potentials up to $ < n 3 ^ -1 . For 
larger concentrations many-chain intermolecular interactions come into play 

The main problem in this approach is how to obtain the intramolecular interactions, 
i.e., the potentials among the blobs of the same chain. The problem is trivial for a 
dumbbell model (n — 2), where the blob-blob interaction is just the potential of mean 
force, but becomes increasingly difficult when increasing the number n of blobs per chain. 
Indeed, the intramolecular effective interaction, which is simply minus the logarithm of 
the joint distribution of the blob positions, is inherently a many-body interaction, which 
cannot be represented as a hierarchy of two-body, three-body, etc. terms, at variance 
with what happens for the intermolecular potentials. In that case the small-density 
expansion gives the pair potential at lowest order, the three-body potential at next-to- 
leading order, and so on. Therefore, approximations must be introduced in trying to 
reproduce some features of the underlying full-monomer system. 

As a first attempt in this direction, Pierleoni et al. [32J introduced a multiblob 
model, referred to as model Ml in the following, for homopolymers in good-solvent 
conditions. They started from the intramolecular and intermolecular potentials 
appropriate for dumbbells (a two-blob molecule), a problem that can be solved as 
explained in Addison et al. [33] The potentials for chains with more blobs were then 
obtained from the dumbbell potentials using a transferability hypothesis. This model 
has the correct scaling behavior for good-solvent polymers in dilute and semidilute 
solutions, [32] including the excluded-volume screening at large length scales with 
density. However, its prediction for the EOS of the solution is incorrect, [31] shading 
serious doubts on the correctness of the transferability assumption for the potentials. A 
modified model, referred to as model M2 in the following, was also introduced. [34] In this 
model, for each number of blobs, the parameters of the potential are tuned to match full- 
monomer results for the EOS. Although successful in reproducing the thermodynamics, 
model M2 is inherently different from Ml in that full-monomer results at finite density 
are needed to tune the model parameters, which is an evident limitation of this approach 
in more complex situations. 

In this paper, beside comparing in detail results for models Ml and M2 to full- 
monomer predictions in a wide range of polymer concentrations in the semidilute regime, 
we introduce a new coarse-grained model for semidilute solutions. In this model we map 
a single chain on a tetramer, i.e., a chain of four blobs, in such a way to reproduce all 
two-body, and some three-body and four-body distributions of an isolated good-solvent 
polymer at the full monomer level. In order to accomplish this program we use bonding, 
bending and torsional angle potentials, plus additional 1-3 and 1-4 central potentials 
which we obtain with the Iterative Boltzmann Inversion (IBI) procedure. [351 EHl 137] 
Furthermore, a single intermolecular pair potential between blobs of different tetramers 
is obtained in such a way to reproduce the center-of-mass radial distribution function 
between two isolated chains at the full-monomer level. This model provides correct 
results for the EOS up to reduced densities $ ~ 2, a considerably larger range of 
densities with respect to simpler models. The tetramer model presented here is a first 
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successful attempt to realize the multiblob program: building a multiblob coarse-grained 
model based on zero-density potentials which is able to provide the correct single-chain 
structure and EOS at finite density. 

Our approach is very close to the coarse-graining procedure applied by Fritz et 
q/.[38] to polystyrene. Similar methods were also used for coarse-grained simulations 
of typical benchmark chains like polycarbonates and polystyrene in a melt, [39] HQ] 
although in these works the potentials were fixed by requiring the coarse-grained model 
to reproduce structural properties at fixed pressure and temperature — hence potentials 
also depend on thermodynamic variables, as is the case for the density- dependent 
potentials. [271 E] In some coarse-graining procedures also thermodynamic information 
was taken into account to fix the potentials, see, e.g., Rossi et a/. [31] and references 
therein. We should also mention the approach of Vettorel et al. [42], which extends 
previous work on a single-blob model. [13] In this multiblob model each blob carries 
internal degrees of freedom, to account for the density profile of the underlying full- 
monomer subchains. However, the interactions are not derived ab initio in the coarse- 
graining procedure, but are obtained by using phenomenological arguments. Finally, 
we mention the work of Clark et al. [H] which applies integral-equation methods to a 
coarse-grained model appropriate for polymer melts. 

The paper is organized as follows. In section [2] we present the general formalism 
behind any coarse-graining procedure and we report our specific methodology to derive 
the tetramer model. In section [3] we compare the structure and the thermodynamic 
behavior predicted by the tetramer model with those (referred to as full-monomer results 
in the following) obtained by using lattice polymer models with a large number of 
monomers (L > 1000) — hence appropriate to obtain the universal, scaling behavior 
- both at zero and at finite density in the semidilute regime. In section [4] we report 
results for the coarse-grained models Ml and M2 and compare them with the tetramer 
and the full-monomer data. Finally, we collect our conclusions and perspectives in the 
last section. In the appendix we give universal predictions for the blob radius of gyration, 
an important quantity to obtain a meaningful comparison between any coarse-grained 
model and the underlying full-monomer model. 

2. The blob model 

In order to obtain the coarse-grained blob model (CGBM), one works in the zero- 
density limit and determines in successive steps the intramolecular potentials, the 
two-body intermolecular potentials, then, at least in principle, the three-body, four- 
body, etc. intermolecular potentials. In an exact mapping all /c-body intermolecular 
interactions should be considered. However, as discussed in the introduction, higher- 
order intermolecular interactions can be neglected if one only considers small densities 
$ ^ $max, where $ max ~ n Zv ~ x increases (for a given level of approximation) with the 
number n of blobs. 
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2.1. The blob representation of the polymer 

In the multiblob approach, the basic object is the "blob", which is a subchain of the 
polymer. Suppose we wish to partition a polymeric chain of L monomers into n blobs 
of m = L/n monomers each. If the monomer positions are given by {ri, . . . , r^}, one 
first defines the blob positions Si, . . . , s n as the centers of mass of the subchains of m 



monomers, i.e. 



j mi 

- E r *- ( 2 ) 



m 

fe=m(i-l)+l 



For the new coarse-grained chain {si, . . . ,s„} one defines several standard quantities. 
First, one defines its radius of gyration 



R 2 



9 > b 2n 2 

ho 

Such a quantity is always smaller than R g , since 

% = K>+ 1 -T,<-i« <«> 

i 

where r 9ti is the radius of gyration of z-th blob: 

mi 

k,l=m(i— 1)+1 

The ratios R gb /R 2 and r 2 g i /R 2 of their averages over the polymer configurations are 
universal, hence independent of the nature of the underlying polymer model as long as 
L is large enough. The average of the blob squared radius of gyration r g defined by 
r 2 = (l/n) X]i r g,i scales quite simply with n in the zero-density limit. As discussed in 
Appendix A[ for all n > 4 we have quite precisely 

^- = 1.06n- 2y , (6) 
R? g 

an expression we will use extensively in the present work (here and in the following we 
will use a hat to indicate zero-density quantities). 

Beside the radius of gyration, we can consider the bond-length distributions (all 
blob distributions depend on the number n of blobs, which is implicit in the notation) 

^(r) = (5(| Sl - Sj |-r)), (7) 

where (•) is the statistical average over all chain conformations, which satisfy the 
normalization conditions 

dr P i3 (r) = 1. (8) 

% Note that here we use the same notation for the squared radius of gyration of a single-chain 
configuration and for its statistical average over all chain conformations. When we will need to 
distinguish between the two quantities, the average squared radius of gyration will be indicated as 

(Rl)- 
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They depend on the chosen length scale. As it is standard in renormalization-group 
analyses of polymer behavior, the relevant quantities are the adimensional combinations 
R g Pij(r). For L — > oo they converge to universal, hence model-independent, functions 
fij(p) with p = r/Rg, which are normalized as 

poo 

/ dpfy(p) = l. (9) 
JO 

Note that, as usual, scaling functions depend only on the adimensional combination 
P = r/R g . 

In this paper we will also consider the adimensional intramolecular distribution 
function 

2R 3 

fetra(r) = — ^- $> (3) (s, - s, - r)>. (10) 



k 3 



For large L, <?mtra(?") converges to a universal function Gmtra(p), P = r /Rg, which is 
related to the bond-distribution functions defined above by 

Note that the the ratio R 2 b /R 2 is simply related to the second moment of gi n t ra ( r ) [ m 
the scaling limit to that of Gintra(p)]- For L — > oo we have 

R lb _(n-l) 



R 2 2n 



p 2 G in Up)d"p- (12) 



Beside two-site distributions, one can define three-site correlation functions 

P M - fc (n,r 2 ) = (5 i3 \s l - sj - r 1 )<j( 3 )(s i - s fc - r 2 )) (13) 

- the corresponding adimensional combinations RgPijk(ri,r 2 ) converge to universal 
functions of Yi/R g and r 2 /R g — and, analogously, four-site, five-site, etc. correlations. 

As a check of the quality of our results we shall often consider the distribution of 
R g> b- More precisely, for each polymer configuration we consider the corresponding 

^1/2 

radius R g>b and the adimensional ratio R gjb /(R 2 ) , where (R 2 ) is the average of 
the squared radius of gyration over the polymer configurations. The corresponding 
distribution 

PrA%) = (s ^-^= - q b j ^ (14) 

is universal in the large-L limit. Note that this distribution function cannot be written 
in terms of the bond-length distributions, but is instead a particular n-blob correlation 
since R g ^ depends on the positions of all blobs. 
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2.2. The coarse-grained model 

In the CGBM the basic object is a polyatomic molecule with n atoms located in 
ti,...,t n . All length scales are expressed in terms of R g , hence potentials and 
distribution functions depend on the adimensional combination p = t/R g . The 
intramolecular potentials are determined by requiring all adimensional distributions 
to be identical in the polymer model and in the CGBM at zero density. For instance, 
we require 

{5(\ti - tjl/Rg - p))ogbm = fM, (15) 

where (-)cgbm is the average over all single-chain CGBM configurations and fij(p) are 
the universal functions defined above, which are computed in the polymer model. 

In principle, the determination of the intramolecular potential is straightforward. 
First, one determines the n-body blob distribution in the polymer model at zero density: 

P n (r 12 , . . . , r ln ) = (f[ 6® (s fc - si - r lfe )\ , (16) 

where the average is over all single-polymer conformations. The adimensional 
combination R g P n converges for L — > oo to a universal distribution: 

^ (n - 1} P n (r 12 , . . . , r ln ) = /„ (p 12 , . . . , Pln ) . (17) 

where p^ = Vij/R g . The CGBM intermolecular potential is then 

0V(pl, ■ ■ ■ , P n ) = ~ lo &fn (P 2 - Pi, ■ ■ ■ > Pn ~ Pi) > ( 18 ) 

where p t = ti/R g . By definition, this choice ensures that the distribution of the n 
atoms in the CGBM is identical to the distribution of the n blobs in the polymer model. 
Hence the intramolecular structure is exactly reproduced. Note that potential f|T8l) is an 
intrinsically n-body interaction and thus there is no natural method to represent it as 
a sum of two-body, three-body, etc., terms. Because of the universality of the function 
f n , the potential is independent of the polymer model and is valid for any polymeric 
system under good-solvent conditions. 

The radius of gyration R 9 ,cgbm of the CGBM molecule differs from the polymer 
radius of gyration R g but agrees instead with R g ^. It is important to take this difference 
into account when comparing finite-density results. For polymers, the behavior is 
universal once densities are expressed in terms of the polymer volume fraction 

where N is the number of polymers contained in the volume V. Full-monomer results 
should be compared with results obtained in the CGBM at volume fractions 

^ = y r *V> m 

where iV& is the number of CGBM molecules. Note that R g and not R g ^ appears in the 
definition of $&. Since R gt b/R g converges to 1 as n increases, for n large, say n > 30, this 
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conceptual difference is not relevant in practice. In our case, instead, since we consider 
n — 4, it is crucial to use the correct definition, that is the quantity 

Once the intramolecular potentials are determined, one must determine the 
intermolecular potentials, which must be such to reproduce the potentials of mean force 
in the polymer model. Note that, in order to have an exact mapping of the polymer 
model onto the CGBM, not only should pair potentials be considered, but also three- 
body, four-body, etc. interactions should be included. [261 EH] However, as we already 
discussed in the introduction, as n increases, these many-body interaction potentials 
are expected to become smaller, so that the CGBM with only pair potentials should be 
accurate in a density interval which widens with increasing n. 



2.3. Determination of the four-blob CGBM intramolecular potentials 

In order to have an exact mapping of the polymeric system onto the n-blob CGBM, 
one should consider an n-body intramolecular potential, which, for n > 2, can be 
expressed in terms of 3(n — 2) scalar combinations of the positions of the blobs because 
of rotational and translational invariance. The complexity increases rapidly with n and 
for this reason we decided to consider the case n = 4, which allows us to limit the 
number of approximations needed and, at the same time, allows us to go beyond the 
dilute regime up to $ ~ 2-3. However, even for n as small as 4, an exact determination of 
the intramolecular potential requires considering a function of 6 independent variables, 
which is far too complex in practice. 

Thus we have used a limited set of interactions. The intramolecular interactions 
have been modelled by introducing six different potentials, each of them depending on 
a single scalar variable. This choice is arbitrary, but, as we will show in the following, 
it is particularly convenient and works quite well. First, we consider a set of bonding 
pair potentials: atoms i and j of the tetramer interact with a pair potential (p) with 
p = |tj — tj\/Rg. Because of symmetry we have Vi 3 (p) = V^p) and Vi 2 (p) = V 3i (p), 
so that there are only four independent potentials to be determined. Then, we consider 
a bending-angle potential Vf,(cos/3) and a torsion-angle potential V t (8), where (3 and 9 
are defined as 

At, ■ At m 

CQS ^ = lAt.nAwr (21) 

= (At, x At m ) ■ (At m x At, +2 ) 
| At, x At m ||At l+1 x At l+2 | 
with Atj = t i+1 — tj. Note that in the tetramer there are two bending angles, which 
are equivalent by symmetry, and a single torsion angle. This particular form of the 
potential set is inspired by the usual modelling of bonded interactions in macromolecules. 
However, in that context one only considers a bonding potential between atoms which 
are first neighbors along the chain, a bending and a torsional term. Instead, our 
parametrization includes interactions between atoms that are not neighbors along the 
chain, thereby taking into account to some extent cross-correlations among different 
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degrees of freedom. We note that the bending potential and the torsion potential involve 
three and four atoms, respectively, and thus allow us to introduce some of the three-body 
and four-body interactions present in the exact parametrization. 

Since we are using a limited set of interactions, not all distributions of the internal 
degrees of freedom can be exactly reproduced by the CGBM. We must therefore choose 
the distributions which we wish to be identical in the polymer and in the CGBM 
case. Given our choice of interaction potentials, it is natural to use the adimensional 
bond- length distributions R g Pij(r) and the distributions of the bending and torsion 
angle [in the blob representation of the polymer model, these angles are defined by 
replacing t with s in Eqs. (121j) and ( 1221) ]. which are particular three-body and four-body 
correlation functions. If we indicate collectively the six potentials to be determined 
with Vi(xi), the (adimensional) distributions of the Xi variables with Pi(xi) in the 
CGBM and with -P^fm^j) in the full- monomer case, the potentials should be such that 
t FM{xi)- The universal (i.e., model-independent) distributions P^FM{xi) in 
the polymer case have been determined by performing simulations of self-avoiding walks 
on a cubic lattice. To detect scaling corrections, we consider chains of length L = 2100, 
4100 (the corresponding blobs have L/4 = 525,1025 monomers, respectively). The 
(adimensional) distributions obtained in the two cases agree within errors, indicating 
the absence of finite-length effects. 

We determine the potentials of the CGBM by using the Iterative Boltzmann 
Inversion (IBI) scheme. [351 EH EZ] m this approach the effective interactions which 
reproduce the target structural quantities are determined iteratively. The potentials of 
mean force of the corresponding full-monomer distribution, — In Pj ^(ij), have been 
chosen as initial guesses for all interactions except for Vu{r). For Vu(r) we have assumed 
a simple Gaussian potential: we use a Gaussian approximation of the potential of mean 
force between two polymer centers of mass, rescaling the width of the Gaussian with 
the ratio of the radii of gyration of the blob and of the entire chain. 

At the end of the optimization procedure, the bond and angle distributions are 
reproduced quite precisely, see figures [TJ and |2j The potentials obtained are reported 
in figure |3j The pair potentials Vi 2 (p) and V^p) are very similar, indicating that 
end effects are not very important. They have a minimum for p m 0.5 and are very 
soft at the origin: V(0) — V(0.5) ~ (0.9-1.0) ksT. For p — \ oo they increase quite 
rapidly and for 2 < p < 3 they approximately behave as p 2 2 (Vi 2 ) and p 2A (V23). 
The potential between next-nearest neighbors is continuously decreasing and hence it 
penalizes configurations in which the two atoms are close. Finally, Vu(p) appears to 
be irrelevant for p > 1. As for the bending- angle potential, it penalizes configurations 
with (3 < 90°, while it is essentially flat for > 90°: the potential has a minimum for 
ps 100° and V b {180°) - V 6 (100°) « 0.Uk B T. Finally, the torsion potential turns out 
to be quite flat: it increases with 6 and changes only by 0.3/csT, going from 0° to 360°. 

The results for the potentials are quite interesting since they indicate the presence 
of an effective hierarchy among the different contributions. The pair potentials between 
neighbors are the most important. For instance, for typical configurations, say, for 
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Figure 1. Adimensional two-body distributions R g Pij(r) as a function of p = r/R g . 
We report full- monomer results (dashed line) and the results for the tetramer CGBM 
with the potentials obtained by means of the IBI procedure. 




Figure 2. Distributions of the torsion angle 9 (right) and of the cosine of the bending 
angle /3 (left) as a function of 9 (in radians) and of cos j3, respectively. We report full- 
monomer results (FM, open symbols and solid line) and the results for the tetramer 
CGBM (t, closed symbols and dashed lines) with the potentials obtained by means of 
the IBI procedure at $ = 0, $ = 1.09, $ = 4.36 and $ = 8.72. For sake of clarity, 
results at different densities are shifted upward by constant values reported in the 
legend. 



0.4 < pi2,P23 ^ 1.5, see figure U the potentials V^ipn) and ^3(^23) vary significantly, 
by 4k B T-5k B T. Instead, for typical distances 0.7 < p± 3 < 2.2, the potential ^13(^13) 
varies much less, approximately by 2/cgT, while in the typical interval pu > 1, Vi^pu) 
varies only by 0.03ksT. Clearly, the relevance of the interactions decreases as the 
chemical distance between the atoms increases, even though interactions between distant 
atoms can never be neglected, otherwise one would model a random-walk chain and not 
a polymer under good-solvent conditions. 

The bending- angle potential varies at most by k B T and is thus less relevant than 
the bonding potentials, while the torsion potential only provides a small correction. 
This seems to indicate that the relevance of the interactions decreases with the number 
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Figure 3. Potentials for the tetramer CGBM. In the left column we show the pair 
potentials between neighbors along the chain (top) and between next-to-nearest atoms 
(bottom), as a function of p = r/R g . In the central column we report the bending-angle 
potential as a function of cos (3 (top) and the torsion-angle potential as a function of 6 
(in radians) (bottom). In the right column we plot the intermolecular potential (top) 
and the potential between the first and last atom of the chain (bottom) as a function 
of p = r/Rg. The symbols represent the numerical results obtained for a discrete set 
of values of p, the lines are interpolations. 



of atoms involved, so that, when the number n of atoms increases, one can safely neglect 
higher-body interactions. 

It is important to stress that the pair potentials ^(p) and V^p) are somewhat 
different from the potentials one would obtain by using the transferability hypothesis 
as suggested by Pierleoni et a/.[32j If we use Eq. ((6]) to relate f g to R g , we would obtain 
a transferability potential (see the expression reported in the caption of figure 1 of 
Pierleoni et al.[32\) 

V tr (p) = 1.92exp(-3.85p 2 ) + 0.534(2. 19 p - 0.73) 2 , (23) 

where p = r/R g . The minimum of this potential occurs for p « 0.67, to be compared 
with p fa 0.5 of V12 and V23. Overlaps are more suppressed since V tr (0) — Vj r (0.67) ~ 
1.55k B T (for our potentials we find 0.9kBT-1.0k B T). Morever, the potential V tr (p) 
increases much less than Vi 2 (p) or V 23 (p) as p increases. For instance, for p 1.5, 
a value which occurs quite frequently, see figure [U we have V tr (p) ~ 3.5k B T, while 
V 12 (p) ^ A.9k B T, V 23 (p) ^ Q.lk B T. 

2.4- Determination of the CGBM intermolecular potentials 

As for the intermolecular potentials, we have made some drastic simplifications. First, 
we do not consider n-body interaction terms, which, as we already mentioned, are 
important only for densities $ > n Zu ~ l . Then, we consider a single intermolecular 
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pair potential W(p): the atoms interact with the same potential, irrespective of their 
positions along the chains. Such a potential has been obtained by requiring the CGBM 
to reproduce the center-of-mass intermolecular distribution function. Indeed, define in 
the polymer model 

9cu{r) = (e-^) 0>1 , (24) 

where (-)o, r indicates the average over two polymers, the centers of mass of which are 
in the origin and in r, respectively, and U12 is the intermolecular energy. In the scaling 
limit L — > 00, gcAi(r) converges to a universal function fcM(p), P = r /Rg- The pair 
potential has been determined so that 

9cm,cgbm(p) = fcAi(p), (25) 

where Qcm<cgbm{ l P) is the corresponding distribution function in the CGBM. Note that 
the second virial coefficient B 2 is related to gcMKj) by 

B2 = \J dh[l - g CM (r)} = 2vr J r 2 dr[l - g CM (r)}. (26) 

Hence, equality (125j) guarantees that the adimensional combination A2 = -B2/-R3, hence 
the thermodynamics in the small-density limit, agrees in the CGBM and in the polymer 
model. 

The potential (3W(p) has been parametrized as 

p\V(p) = Cl exp(-c 2 p 2 ), (27) 

in terms of two unknown parameters c± and c^. They have been determined following 
the approach of Akkermans et a/.[IHJ 06] Requiring the model to reproduce at best 
the polymer scaling function p 2 fcAi(p), we obtain the optimal values c\ = 1.66 and 
C2 = 3.9. For these parameter values the model with potential (12 7p has an intermolecular 
pair distribution function which agrees quite precisely with the corresponding polymer 
quantity, see figure HI The result depends on the parametrization and we cannot exclude 
that a different parametrization with the same number of parameters gives results of 
better quality. Potential (1271) differs only slightly from the intramolecular potential 
Vw{p) — 1.86 exp(— 4.08425p 2 ). Interactions between the atoms at the ends of the chain 
or between atoms that belong to different chains are quite similar. It is interesting to 
compare our result with that one would obtain by using the transferability hypothesis 
as suggested by Pierleoni et al: [32J (3W tr (p) = 1.92 exp(— 3.85p 2 ). The range of the 
potential is the same, but the potential we obtain is somewhat softer. 



3. Comparison of CGBM and polymer results 

In order to understand how well the tetramer model reproduces the polymer behavior 
we have performed extensive simulations of the tetramer and of a polymer model at 
zero density and at volume fractions $ = 1.09, 4.36, 8.72. Since CGBM and polymer 
results should be compared taking <I> = <E> 6 , see (fT9|) and (!20l) . we drop the suffix and 
thus <l> also refers to At finite density polymers have been modelled by means of the 
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Figure 4. Correlation function gcAi{p) as a function of p = r/R g . We report full- 
monomer results (dashed line) and results (filled circles) for the model with potential 



Domb- Joyce (DJ) model with w = 0.505838, a value which is close to the optimal one 
for which no leading scaling corrections are present (see Caracciolo et al.\A7\ for details 
on the model). It allows us to determine precisely the universal, model- independent 
scaling functions by using chains of moderate length. We consider walks of length 
L = 600 and L = 2400, verifying the practical absence of scaling corrections. As in 
previous work [48], we considered different box sizes, finding negligible size effects when 
the number of chains in the box exceeds 100. Simulations have been performed using 
the algorithm described in Pelissetto. jlH] 

3.1. Zero-density 

By construction, the tetramer CGBM reproduces the bond-length distributions. As 
we have already remarked in Sec. 12.11 see (fill and (112"]) . the ratio R g ^/R g can be 
expressed in terms of these distributions, hence this ratio should assume the same value 
in the tetramer and in the polymer case. Numerically, we find R g ^/R g = 0.89093(7) 
for the tetramer and R g) b/R g = 0.89210(11) for the polymer case. The difference 
is approximately 0.1%, which shows how accurate the intramolecular potentials we 
determined are. Moreover, not only R g>b /R g agrees, but also its distribution f }14"|) is 
the same for polymers and tetramers, see figured). This is a nontrivial check, since 
this distribution is not directly related to the bond-length distributions, nor to those of 
the bending and torsion angles. Clearly, the tetramer models correctly the shape and 
size of a polymer at zero density. 

Since we have matched the center-of-mass distribution function to determine the 
intermolecular potential, the tetramer CGBM should give the correct second virial 
coefficient. If we expand the compressibility factor as 

Z = y^r = 1 + B 2 c + B 3 c 2 + 0(c 3 ), (28) 




Figure 5. Adimensional distribution PRfi{%) of the ratio <j<& = Rg^/iRg) 1 ^ 2 , see 
definition (Tl4|) . for the tetramer (t) and for polymers (FM). We report full-monomer 
(FM) and tetramer (t) results for $ = (a), $ = 1.09 (b) and $ = 8.72 (c). 

the quantity A 2 = B 2 /R? g is universal. An accurate estimate is jU] A 2 = 5.500(3). 
For the tetramer we obtain A 2tt = 5.597(1). The difference is approximatey 1.8% and 
is representative of the level of precision with which the tetramer model reproduces 
the center-of-mass distribution function. Much more interesting is the comparison of 
the third virial coefficient, since it provides an indication of the accuracy with which 
the tetramer model reproduces the polymer thermodynamics in the dilute regime and 
also of the importance of the neglected three-body forces. The universal combination 
A 3 = B 3 /R® was computed by Caracciolo et al. [4T] finding A 3 = 9.80(2). 

In order to determine A 3 , two contributions had to be computed. One contribution 
is the standard one, the only present in monoatomic fluids and in fluids of rigid 
molecules, A' 3 « 10.64, while the second one is a flexibility contribution A 3 ji « —0.84 
(it corresponds to —T\R~ & in the notations of Caracciolo et a/. [47]). The combination 
A 3 as well as the two contributions A' 3 and A 3 ji are universal, hence it makes sense to 
compare them with the tetramer corresponding results. We obtain 

A 3)t = 9.99(2), A' 3it = 10.57(2), A w = -0.581(5). (29) 

The tetramer reproduces quite reasonably the third virial coefficient, the difference 
being approximately 2%. Note that much of the discrepancy is due to A 3 jf. the 
tetramer is more rigid than the polymer. It is interesting to compare these results 
with those obtained by using the single-blob model in which polymers are represented 
by monoatomic molecules interacting by means of density-independent potentials, [jj 
If we use the accurate pair potential of Pelissetto et a/.[llj we obtain A 3 = 7.844(6) 
(of course here A 3 ji = 0) which deviates by 20% from the polymer result. Hence, the 
tetramer model represents a significant improvement for the thermodynamics. 

§ If we were using density-dependent potentials, the thermodynamics, hence all virial coefficients, would 
be exactly reproduced. However, here we only consider models with density-independent potentials, 
hence the only meaningful comparison is with the single-blob model in which the interactions are density 
independent. 
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Figure 6. Three-body potential of mean force /3V3(ri2, ri3, r23) for ru = ri3 = ^23 = 
r, as a function of p — r/R g . We report the tetramer (t) result and the prediction of 
full-monomer (FM) simulations. 



As a further check we compute the effective three-body potential of mean force 
defined by [261 H] 

Fl3 ' F23) = " M (e-^") ri2 (e-^) ri3 (eX )r2 ; ( 30 ) 
here Uij is the intermolecular potential energy between tetramers i and j and the average 
(')ri2,ri3,r 2 3 is over triplets of tetramers such that = — r^, where rj is the position 
of the center of mass of tetramer i. 

We computed (3V 3 (r 12 , r 13 , r 2 3) for equilateral triangular configurations such that 
r i2 — r i3 — r 23 — r - The result is reported in figure [6] and compared with the analogous 
quantity computed in full-monomer simulations. At variance with the single-blob model 
for which f3V% = 0, the tetramer model reproduces the polymer j3V 3 quite reasonably: 
differences — the tetramer potential is slightly more attractive - - are observed for 
p = r/R g < 1, but they are only significant for p < 0.5, i.e., when the tetramers are 
very close. This is consistent with the analysis of the third virial coefficient: in the dilute 
limit three-body interactions are correctly reproduced by the tetramer model, without 
the need of introducing a three-body potential among tetramer atoms. 



3.2. Semidilute regime 

As we have discussed in the introduction, the tetramer model is expected to reproduce 
the full- monomer results up to $ ~ 2, representing a significant improvement with 
respect to the single-blob model which shows large deviations already for $ = 1. To 
check this expected behavior we compare tetramer and full-monomer simulations at 
$ = 1.09, 4.36, and 8.72. 

Let us begin with the structural properties. In figure [7] we report the adimensional 
intramolecular distribution function g\-at Ta ,{r) . For $ = 1.09 the agreement between the 
tetramer and the full monomer results is excellent. However, as $ increases, deviations 
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Figure 7. Adimensional intramolecular distribution function (/intra (j) as a function 
of p = r/Rg. We report tetramer (dashed lines, t) and full-monomer (full lines, FM) 
results. Distributions corresponding to $ = 4.36 and $ = 8.72 are shifted upward for 
clarity according to the legend. 
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are observed for p — r/R g < 1. For $ > 4 the tetramer is more swollen than the polymer: 
the probability for two blobs to be at a given distance p < 1 is significantly smaller in 
the tetramer than in the full-monomer chain. These results are further confirmed by 
the results for the radius of gyration. For the tetramer we have 

0.85636(4) $ = 1.09, 

0.8181(2) $ = 4.36, (31) 

0.8047(1) $ = 8.72, 
to be compared with the full-monomer results 

0.8523(2) $ = 1.09, 

0.7823(2) $ = 4.36, (32) 

0.7346(6) $ = 8.72. 

For $ = 1.09 the agreement is very good, consistently with the results reported in 
figure As $ increases, however, the tetramer is more rigid than the polymer and 
Rg,b{&) I ' Rg is larger in the tetramer than in the polymer case. The same conclusions 
are reached by looking at the distribution of R 9t b, see figure El For $ = 1.09 the 
agreement is excellent, while for $ = 8.72 the tetramer distribution is slightly shifted 
towards larger values of R g ^- 

It is also interesting to compare the results for the bending and torsion angles 
reported in figure EJ The distributions appear to have a tiny dependence on $ and to be 
reasonably reproduced by the tetramer for all values of $. For instance, for the largest 
value of $, $ = 8.72, we have for polymers P b (cosf3 = -1) « 0.93, P t (6 = 0) « 0.346, 
to be compared with 0.88 and 0.354, respectively, in the tetramer case. 

Let us now consider the thermodynamics. For this purpose we computed the 
compressibility factor Z = (311/ c for the tetramer using the molecular virial method 



Coarse-graining strategies in polymer solutions 



18 



Table 1. Compressibility factor for the tetramer model, Z t ($), and for polymers, 
Z p ($). Polymer results are taken from Pelissetto [48]. 



$ 


z t ($) 




0.054 


1.07363(3) 


1.0725 


0.135 


1.18993(4) 


1.1871 


0.27 


1.39852(6) 


1.3929 


0.54 


1.8499(1) 


1.8536 


1.09 


2.9090(1) 


2.9589 


2.18 


5.2660(2) 


5.6342 


4.35 


10.2056(4) 


12.229 


6.53 


15.2279(1) 


20.019 


8.72 


20.2811(1) 


28.716 




Figure 8. Compressibility factor for polymers (FM), for the tetramer (t) and for 
the single-blob (mono) coarse-grained model. For the tetramer we also show the 
interpolation Z($) = (1 + + a 2 <S> 2 + a 3 $ 3 ) 1 / 2 /(l + a 4 $) 1 / 2 , with ai = 3.32676, 
a 2 = 4.67289, a 3 = 3.58551, a 4 = 0.65439, of the data reported in Table HJ 

[4"9[ [50] (c = N/V is the number concentration). As Z is dimensionless, polymer and 
tetramer results at the same value of $ can be directly compared. Estimates are reported 
in Table [TJ For $ < 1 the tetramer Z is very close to the polymer prediction: for 
$ = 1.09 it differs by 2% from the correct result. As $ increases, however, differences 
increase and the tetramer model underestimates the correct pressure. In figure M we 
compare Z(<p) for the tetramer with the corresponding expression for polymers. |48j At 
the scale of the figure, good agreement is observed up to $ ~ 2. For larger densities, 
Z(<$>) in the tetramer increases slower than in the polymer case. Indeed, while in 
polymers we expect Z ~ (frVC 3 ^ 1 ) ~ for large for the tetramer Z is expected to 
increase only linearly with $ (since the potential is soft, for $ — >■ oo the random-phase 
approximation should become exact [5T]). 
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Figure 9. Center-of-mass distribution function gcAi(p) as a function of p = r/R g 
for polymers (FM, open symbols) (obtained from full- monomer simulations), for the 
single-blob coarse-grained model (Mono, left, closed symbols) and for the tetramer (t, 
right, closed symbols). We report results for <£> = 1.09 and $ = 4.36. 



As can be seen from figure El the tetramer model is significantly better than the 
single-blob modeljJJ] in which each polymer is represented by a single atom. At $ = 1.09 
such a model gives Z = 2.70(1), which underestimates Z by 8%, much more than the 
tetramer model (at this density the error on Z is 2%). The single-blob model gives a 
value for Z which differs from the polymer one by less than 2% only up to $ ~ 0.38, 
i.e. up to densities which are a factor-of-three smaller than the corrisponding ones for 
the tetramer. This improvement confirms the scaling argument we presented in the 
introduction. As we explained in the introduction, the multiblob model should give the 
correct thermodynamics up to a polymer volume fraction $ max which scales as n? v ~ x . 
If we compare the tetramer model with the single-blob one, we thus expect the density 
range in which the model is predictive to increase by 4 3iy ™ 1 « 2.9, which is indeed what 
we find. 

Let us finally compare the center-of-mass distribution function gcNi(p)- It is 
reported in figure for $ = 1.09 and $ = 4.36. In the first case, the tetramer result 
is on top of the polymer results. For $ = 4.36 small discrepancies at short distance 
(p ~ 0-5) are present. For instance, for the tetramer we have <7cm(0) = 0.591(6), 
<?cm(0.5) = 0.688(1) at $ = 4.36, respectively, to be compared with the polymer results 
<7cm(0) = 0.550(5), gcA/(0.5) = 0.660(2). These small differences are responsible for the 

| For the single-blob model we shall always use the accurate expression of the pair potential given in 
Pelissetto et al. [TI] 
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differences in the estimates of Z since gcu is related to Z by the compressibility rule^jl 

^) =l + cj(gc M (r)-l)dh. (33) 

Note that, even though the thermodynamics is poorly reproduced, also the single-blob 
model gives an estimate of <?cm ( r ) which is only slightly different from the polymer one. 
The largest differences are observed for p — > 0. At overlap we obtain <7cm(0) = 0.344(1) 
and (?cAf(0) = 0.613(1) for $ = 1.09 and $ = 4.36, to be compared with the polymer 
results 0.307(9) and 0.550(5). 



4. Comparison with other models 

In the previous section, we discussed the finite-density behavior of the tetramer model 
and we showed that it is quite accurate, for both structural and thermodynamic 
properties, up to $ w 2, in agreement with the multiblob argument of Pierleoni et 
al. [32J presented in the introduction. It represents a significant improvement with 
respect to the single-blob model, which is unable to reproduce structural properties 
and reproduces the thermodynamics only deep in the dilute regime (the compressibility 
factor Z for the single-blob model differs from the polymer one by less than 5% only for 
$<0.75). 

Here we wish to investigate the structural and thermodynamic behavior of two 
other models discussed in the literature. 



4-1. Definition of the models 

First, we consider the model introduced by Pierleoni et al. [32J - we name it model 
Ml. A CGBM with n blobs is a chain in which neighboring atoms belonging to the 
same chain interact with an intramolecular potential 

VWd(r) = Ae~ ar2 ^ + k(r/r g - r ) 2 ; (34) 

atoms that belong to the same chain but are not neighbors, or belong to different chains 
interact with a potential given by 

VU-bond(r) = Ae- ar2 " 2 ^ (35) 

where A and a are the same as in fl34|) . In these expressions f g is the average zero-density 
radius of gyration of the blobs and sets the length scale. The model depends on several 
constants, which can be easily determined in the dimer case, i.e., for n = 2 (see caption 

% In principle one can use this expression to determine the density derivative of cZ. However, since the 
relevant length scale for gcAi(p) is the average distance d between the centers of mass of the polymer 
chains, finite-size effects will be small only if d/M -C 1, where M is the size of the box containing the 
system. On the other hand, if one uses the intermolecular structure factor, the relevant scale is the 
correlation length £, which, in the semidilute regime, is significantly smaller than d. Therefore, £/M 
is smaller than d/M, which implies that determinations using the intermolecular structure factor show 
smaller finite-size effects than those using <?cm(p)- 
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of figure 1 in Pierleoni et a/. [32]): A = 1.92, a = 0.80, k = 0.534, and r = 0.730. 
To extend the model to values n > 2, Pierleoni et al. [32] made the transferability 
hypothesis: equations ( 1341) and ( 1351) hold for any n, the n-dependence being completely 
taken into account by the radius of gyration of the blob. 

As discussed in Sec. 12.21 when comparing the CGBM results to the polymer ones, 
one should use the radius of gyration R g of the reference polymer model and not the 
radius of gyration R gb of the CGBM. The radius R g (or rather the ratio R g /r g , since 
f g is the basic length scale in this approach) can be determined by using two different 
routes. As suggested by Pierleoni et al, [32] one can determine R g ^/f g for the CGBM 
and then use (T4|). Alternatively one can use fl6]). If the model were a good CGBM, these 
two definitions would give the same result, and indeed in the tetramer case they do. 
Instead, for model Ml we observe quite large differences. For instance, for n = 30, we 
find R 2 ft 41.8fJ if we use R 2 = R gb + f 2 and R 2 g ft 51Af 2 , if we use R 2 g /r 2 = n 2u /1.06. 
These differences do not disappear as n — > oo. An analysis of Ml results with n < 600 
gives for n — > oo the scaling behavior 

R 2 

= An 2u , A = 0.78(3), (36) 

r g 

which is not compatible with (jBJ). 

In this paper we compare the results obtained by using three different "polymer" 



radii of gyration: 

Rg,i/f g = R g ,b/r g , (37) 

R 9 ,2/f g = y^y^TT, (38) 

R 9 ,s/r 9 = n"/VT06. (39) 



Note that, for large n, definitions R g< i and R 9) 2 are equivalent. On the other hand, as 
we discussed, definition R g ^ differs significantly from the others for any n, including 
the limit n — > oo. Recently, Coluzza et al. [22] suggested that model Ml should 
not be considered as a CGBM, but rather as a generic polymer model in good-solvent 
conditions, so that R g ^ should be used as reference scale. In the following we shall mainly 
focus on the first and third definition and we shall label the corresponding results as 
(Ml,l) and (Ml,3), respectively. 

A proper definition of R g is relevant for two purposes: first, length distributions 
are universal only if one expresses the lengths in terms of R g (the relevant scale is 
p = r/Rg); second, at finite density results should be compared at the same value of 
the polymer volume fraction $ defined in fTT9l . Changing the definition of R g changes 
the definition of both p and hence it is crucial to specify which R g one is using. 
Note that Coluzza et al. [22] introduced a complex procedure to compare CGBM and 
polymer results. Their procedure is fully equivalent to the one we have discussed above: 
to analyze finite-density results, one must compare the results at the same value of the 
adimensional volume fraction $. 
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Figure 10. Intramolecular distribution function p 2 gi n tra(^) versus p — r/R g for 
<J) = 1.09 (left) and $ = 8.72 (right). We report full-monomer results (FM) and 
results for models Ml and M2 with n = 30. In models Ml and M2 (but not in the 
polymer case) we use R g ,b as radius of gyration, both in the definition of p and in that 
of 

The thermodynamic behavior of model Ml was studied by Pelissetto. [34J If R g ^ 
is used as reference scale as recently suggested by Coluzza et ai, [52] the model fails 
to reproduce the thermodynamics unless n is larger than 10 3 , but of course, for such 
values of n, there are several other models — the lattice Domb- Joyce model we use 
is one of them — which better reproduce the universal polymer behavior both for the 
thermodynamics and the structural properties. For instance, for n = 100, which is a 
relatively large value, model Ml overestimates the second virial coefficient combination 
A 2 by 20%. A more detailed discussion will be presented below. 

We shall also consider a second coarse-grained model, [31] we call it model M2. 
Conceptually, this was not intended to be a CGBM, but rather a polymeric model 
which reproduces the asymptotic (number of monomers n — > oo) behavior already for 
small values of n. The n-dependent potentials were tuned so that thermodynamics was 
reproduced for <3> < 10. For n = 26, thermodynamics was reproduced taking potentials 
of the form (I34j) and (135 j) with (f g is no longer the blob size, but simply sets the length 
scale) A = 8.28, a — 1, k — 0.15 and r = 0.653. It is important to stress that R g ^ 
was used as reference length in the optimization procedure employed to determine the 
optimal parameters. Therefore, for consistency, for this model it makes no sense to use 
definitions R g ^ and R g ,3- Hence, whenever we consider model M2, R g should always be 
identified with R g ^. 

We have performed simulations for model Ml for n — 8, 16, 30, 60 and of model 
M2 for n = 30. In the second case, one should in principle compute the appropriate 
parameters for n = 30. We will use here the coefficients computed for n = 26, as we 
expect the changes necessary as n goes from 26 to 30 to be tiny. 

4-2. Numerical results and discussion 

Let us first discuss the structural properties, considering the intramolecular distribution 
function gi ntra ,(r). If we consider models Ml and M2 as generic polymer models — in this 




Figure 11. Intramolecular distribution function p 2 gi n tra(^) versus p = r/R g for 
$ = 1.09 (left) and $ = 8.72 (right). We report full-monomer results for a 30-blob 
representation of the polymer (FM) and results for model Ml with n — 30. In model 
Ml we use R g .3, see definition (|39"|) , as radius of gyration. 
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Figure 12. Adimensional distribution Pr^) of the ratio qt- We report: full-monomer 
results for qb = Rg^R 2 ) 1 / 2 (FM); results for model M2 with n — 30, defining <I> using 
Rgfi, with (ft = Rg.b/ (R 2 b) 1 ^ 2 ; results for model Ml with n = 30. In this last case we 
show two different quantities: in case (Ml,l) we use the same definitions as for model 
M2, while in case (Ml, 3), we set qb = Rg,b/ {Ry^) 1 ' 2 an d define $ in terms of R g ,3- In 
the inset we report again the FM data together with the distribution of qb defined in 
(Tl4]) . as appropriate for a 30-blob system (FM,b). 



case we should use R g ^ as length scale — the corresponding results should be compared 
with the monomer intramolecular distribution function, which is defined in ( ITUl) . taking 
n = L. Estimates of p 2 <?mtra(p), which represents the average distribution of the bond 
lengths, are shown in figure ITUl At $ = 1.09, we observe a reasonable agreement for 
both models: they appear to be able to reproduce the structural properties in the dilute 
regime. For $ = 8.72 model Ml shows some, but still reasonably small, differences for 
0-4 ;$ P % 1- Model M2 appears to be slightly in better agreement, except for small 
P<0.2. 

For model Ml, we can also consider R g ,3, see definition (1391) . for the zero-density 
radius of gyration. In doing this, we implicitly assume that model Ml is a CGBM, as 
the tetramer model, and not just a generic good-solvent polymer model. In figure (TT] 
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Figure 13. Bending (left) and torsion (right) angle distributions. We report full- 
monomer results for a 30-blob representation of the polymer (FM) and results for 
models Ml and M2 with n = 30. For Ml and M2, the polymer volume fraction $ is 
computed by using R g .b- Distributions for $ = 8.72 are shifted for clarity. 



we report the corresponding adimensional intramolecular distribution function, which 
should be compared in this case with the polymer results for a 30-blob coarse-grained 
representation (data labelled FM). Discrepancies are significantly larger than in figure 
ITUl Clearly, structural properties are much better reproduced if R g ^ is used as radius of 
gyration, in agreement with the conclusions of Coluzza et al. [52] Note that if one uses 
R g> b the Ml distributions agree exactly with the polymer ones for n — > oo. Indeed, model 
Ml is, in this limit, a generic polymer model and all models have the same infinite-length 
behavior as long as the same adimensional quantities are compared. As a consequence, 
the discrepancies we observe in figure [TT] do not decrease as n increases. Similar 
conclusions are reached by considering the distribution of the radius of gyration, see 
figure [12j Depending on the interpretation of the models as generic polymer models or 
as CGBMs, one should compare the results with the polymer distributions of R g j (Rg) 1 ^ 2 
or of R 9t b/ (Rg) 1 ^ 2 , where R g ^ is the radius of gyration of a 30-blob representation of the 
polymer chain. However, as shown in the inset of figure [121 the two distributions are 
identical on the scale of the figure, so that this conceptual difference is not relevant in 
practice. Model M2 appears to be the one which gives the best agreement, but, if R 9i b 
is used as a reference scale, also the model-Mi distribution is close to the full-monomer 
one. If instead R g>s is used for model Ml, discrepancies are quite large. 

If model Ml is considered as a CGBM, it makes also sense to compare the 
bending and torsion angle distributions. The results, reported in figure [131 (similar 
to those observed in model M2) have little relation with what is observed for the 
polymer case (the full-monomer distributions we report are those appropriate for a 
30-blob representation of the polymer). Hence, even if the bond- length distributions are 
approximately reproduced, correlations between different bonds, for instance angular 
distributions, are not, and the true polymer shape is quite different from that predicted 
by model Ml. 

Let us now come to the thermodynamics. For both models we have determined 
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Table 2. Estimates of A2 = -E>2-Rg and A3 = B^R g using the different definitions 
of R g for model Ml. Numerically, we find R^Jf" 2 , = 7.987(3), 18.82(1), 40.83(4), 
95.37(3), for n = 8, 16, 30, 60. The universal asymptotic values for polymers are [57] 
A 2 = 5.500(3), A 3 = 9.80(2). 



n 




M 






M 






Rg,b 


R g,2 


^9,3 


Rg,b 


R g,2 


R g,3 


8 


9.225(7) 


7.729(5) 


5.815(1) 


32.0(5) 


22.4(3) 


12.7(1) 


16 


8.258(9) 


7.640(8) 


5.548(1) 


26.0(5) 


22.2(4) 


11.7(1) 


30 


7.55(1) 


7.28(1) 


5.354(2) 


21.4(7) 


20.0(7) 


10.8(2) 


60 


6.95(1) 


6.84(1) 


5.183(6) 


17.9(8) 


17.3(7) 


10.0(4) 




Figure 14. Center-of-mass pair distribution function gcM(r) as a function of 
p = r/Rg for polymers (FM), tetramers (t), and model Ml for n = 8, 16,30 at zero 
density. For model Ml, R 2 g = R 2 g 3 = n 2v /1.06 in the main panel, while in the inset 
we report results (n = 30 only) with R g — R g j, (Ml,l) and R g — R g ^ [see definition 
(155)) ] (Ml, 2). For polymers and tetramers, R g is always the radius of gyration. 

the second virial coefficient B2 and the adimensional combination A2 = I?2-R<^ 3 - The 
parameters of model M2 were determined in such a way to reproduce A2 = 5.500, the 
correct result for infinitely long polymers, hence M2 gives the correct thermodynamics in 
the zero-density limit. The results for model Ml are reported in Table [2] for each choice 
of R g . As already discussed by Pelissetto, [M] if R g ^ is used, A 2 differs significantly 
from the asymptotic result, even for n = 60. If R 9j2 is used, discrepancies are smaller 
for n — 8, but substantially the same for n > 30 (not surprising, since R g ^/ Rg,b — * 1 
as n — > 00). Definition i? g 3 gives apparently better results, but we believe that this 
apparent agreement is fortuitous. Indeed, as n increases, B 2 R g l should monotonically 
decrease, increasing the discrepancy with the polymer case. It is easy to compute the 
asymptotic value. For large n model Ml is a generic good-solvent polymer model, hence 
standard universality arguments predict that B 2 R~\ should converge to 5.500, the result 
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Figure 15. Three-body potential of mean force [3Vz{r, r, r) for equilateral triangular 
configurations as a function of p = r/R g for polymers (FM), tetramers (t), and model 
Ml for n = 8, 16, 30 at zero density. For model Ml, R g = R g ^ in the main panel, while 
in the inset we report results with R g = R g ,b (Ml,l) and R g = R g ,2 [see definition 
(Ml, 2). For polymers and tetramers R g is always the radius of gyration. 
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Figure 16. Left: compressibility factor Z(Q>) for polymers (FM) and model Ml: in 
the latter case $ is defined by using R g , b (data labelled Ml.l) or by using R g ^ (data 
labelled Ml, 3). On the right we show the deviations of the (Ml, 3) results from the 
polymer ones, AZ = 100[Zmi,3($)/£fm($) - !]• 



obtained for infinitely long polymers. |47J. Using ( 1361) we obtain for n — > oo 

B 2 R~l = (l.Q6A) 3/2 B 2 R~* = 5.500(1.06 x 0.78) 3/2 w 4.13, (40) 

which differs by 25% from the correct result. To further confirm that there is 
nothing fundamental in the observed agreement, we plot the zero-density center-of-mass 
distribution function gcNi{p) in figure HH For all values of n it differs significantly from 
the polymer one. In particular, the correlation hole <?cm(0), which does not depend on 
the choice of R g , is significantly deeper in model Ml than for good-solvent polymers in 
the scaling limit. 

In Table [2] we also report the third- virial combination A%. If R g ^ is used, results 
differ roughly by a factor of two from the polymer ones. Discrepancies decrease if R g 3 
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Figure 17. Center-of-mass distribution function gcM{r) as a function of p = r/R g 
for polymers (FM), model Ml (left) and model M2 (right), at densities $ = 1.09 and 
4.36. The number of blobs is n = 30. For models Ml (case Ml,l) and M2 (but not in 
the polymer case) we use R g ^ as radius of gyration, both in the definition of p and in 
that of $. Data labelled (Ml, 3) are obtained by using R g ^ as radius of gyration, both 
in the definition of p and in that of $. 



is used, but again this is accidental. The same argument given above for A 2 shows 
that the combination BsR~l converges to 5.5 for n — > oo, which is roughly a factor- 
of-two smaller than the correct result [47J A 3 = 9.80(2). We also report, see figure \15\ 
the three-body potential of mean force for three chains on an equilateral triangle. We 
observe significant discrepancies: results are significantly worse than those obtained by 
using the tetramer CGBM. 

Let us now compare the thermodynamics at finite density. In figure IT6l we compare 
the compressibility factor for polymers (data labelled by FM taken from Pelissetto [48J) 
with that for model Ml. As expected on the basis of the zero-density results, if R g ^ 
is used in the definition of $, very large discrepancies are observed. Moreover, also 
the dependence on $ is incorrect: Z(<&) increases as $ L13 for 6 < p < 9, which differs 
significantly from the correct scaling $ L31 . Discrepancies are significantly smaller (12% 
at most, see the right panel of figure [16]) if R g ^ is used. Again the agreement appears 
to be accidental, since the center-of-mass distribution function differs significantly from 
the polymer one, see figure [T71 Even worse, for p < 1, results using R gh appear to be 
closer to the correct full- monomer results than those obtained by using R g $. Again, note 
that only if R gb is used the distribution function gcM(p) computed in model Ml will 
converge to the full- monomer one for n — > oo. If R g 3 is used instead, the correlation hole 
is always (even for n — > oo) deeper for model Ml than for true polymers at any given 
value of the polymer packing fraction $ ^ 0. Moreover, ^cm(p) shows more curvature, 
reaching approximately 1 at a slightly smaller value of p. 

By construction, model M2 reproduces the thermodynamics up to $ = 10. Indeed, 
the parameters were fixed by requiring Z(<&) to be equal to the polymer compressibility 
in the dilute limit and for $ = 10. Note that it also gives the correct intermolecular 
pair distribution function, see figure [T7J a result which is not a priori obvious. 
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5. Conclusions 

In the last two decades (but the first proposals [53] can be traced back to the '50s) 
several coarse-grained models have been proposed for polymers in solution under good- 
or ^-solvent conditions. In the simplest approaches polymer chains are mapped onto 
single atoms interacting by means of soft potentials. These classes of models are however 
unable to reproduce the structural properties and give the correct thermodynamics only 
in the dilute limit. To go to higher densities, density-dependent potentials [27J E] may 
be used. However, their determination requires in any case finite-density full-monomer 
simulations, which is what one would like to avoid by using coarse-grained models. 
Moreover, it is not clear how accurate they are in more complex situations in which 
polymers only constitute one species in the solution. To overcome these difficulties, the 
multiblob approach was recently proposed, [32] in which each polymer chain is mapped 
onto a short linear chain of n blobs. This model retains some degrees of freedom and 
thus it should allow us to obtain the correct thermodynamics even in the semidilute 
regime. The main difficulty of this approach is the derivation of the intramolecular 
interactions. In Pierleoni et al. [32] potentials were obtained for any value of n on the 
basis of a transferability hypothesis. However, later [31] it was shown that the resulting 
model did not have the correct thermodynamic behavior, indicating that much more 
work was needed to determine the intramolecular interactions. 

In this paper we consider again the multiblob approach, determining the 
intramolecular interactions by matching universal zero-density polymer distributions^ 
We map polymer coils onto four-blob chains (tetramers) which interact be means of 
bonding, bending and torsional angle potentials. They are obtained by requiring the 
bond-length distributions and the distributions of the bending and torsion angles to be 
the same in the tetramer and in the full-monomer model at zero density. As for the 
intermolecular interactions, we only consider pairwise blob-blob interactions, neglecting 
many-blob potentials. This limits the validity of the model to the regime in which 
blob-blob overlaps are rare, i.e., to blob volume fractions rjb = Cb/cl < 1 [q, is the blob 
concentration and c£ = 3/(47rf|)]. For the tetramer this gives $ < n 3l/_1 w 2.9. 

The tetramer model turns out to be quite accurate up to $ ~ 2, in agreement with 
the argument given above. In this range of densities structural properties as well as the 
thermodynamics are correctly reproduced. For instance, for $ = 2.18 the error on Z(f&) 
is 7%. If we compare the compressibility factor computed in the tetramer model to 
that determined in the single-blob model we observe a factor-of-three improvement, 
indicating that the ideas behind the multiblob approach really work. For $ > 2 
significant deviations are observed, both for the structure — tetramers are too rigid 
- and for the thermodynamics — in the tetramer model becomes significanly 

smaller than for polymers as $ increases. 

+ The polymer distributions are computed by using a lattice model. However, standard 
renormalization-group arguments allow us to conclude that exactly the same results would be obtained 
in the limit L — ¥ oo by using any other — discrete or continuous — model. 
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We have investigated again the model proposed in [52], model Ml, studying in 
detail structure and thermodynamics. We find that the model cannot be considered 
as a consistent CGBM, but should rather be thought as a generic polymer model, as 
recently suggested by Coluzza et al. [52] Since the thermodynamics is poorly reproduced 
for small values of n (we mainly investigate the case n = 30), it is not a surprise 
that for these numbers of blobs intermolecular correlations are significantly different 
from those determined in full-monomer simulations with a large number of monomers. 
On the other hand, internal bond distributions are quite well reproduced. Clearly, for 
small values of n, in spite of the poor thermodynamic behavior, model Ml is able to 
model correctly some features of the polymer shape, though not all of them - - for 
instance, angle distributions are not reproduced. This is consistent with the results of 
Coluzza et al. [52] They studied the geometric structure of polymer brushes, comparing 
results obtained in full-monomer simulations and in model Ml. Also in that case, good 
agreement was observed for some structural properties. 

Finally, we consider the model proposed by Pelissetto.[34j In this case, parameters 
were tuned so that the thermodynamics was exactly reproduced up to $ = 10. We find 
that it also reproduces well the intermolecular structure: the polymer center-of-mass 
distribution function is correctly reproduced in the whole density range $ < 10. As for 
the intramolecular structure, we find that the model gives results analogous to those 
obtained for model Ml. Bond-length distributions are approximately reproduced in the 
density range we have investigated, indicating that also this model correctly reproduces 
some features of the polymer shape. 

In conclusion, we have shown that the newly proposed tetramer model is a 
significant step forward in the development of a consistent coarse-grained model of 
polymer chains based on zero-density interactions. To investigate the semidilute regime 
for large densities, i.e., for $ > 2, multiblob models with n > 4 must be developed. 
In this respect, the most important lesson of the present work is that many-body 
intramolecular interactions cannot be completely neglected, if one aims at a consistent 
multiblob model; their absence in model Ml is probably the cause of its failure in 
reproducing the thermodynamics of polymer solutions. Finally, it would be very 
important — we leave it for future work — to develop an analogous coarse-graining 
strategy for chains in 9 conditions. Here single-blob models with pairwise intermolecular 
interactions fail since thermodynamic stability is only obtained by taking into account 
three-chain interactions. Since the tetramer model reproduces quite nicely three-chain 
correlations in the good-solvent regime, it is the good candidate to attack this problem. 
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Figure Al. Combination n u Qi ~ n v r g: i/R g as a function of (i — l/2)/n for 
n = 4, 10, 20, 30. On the left we report the results for $ = 0, on the right for $ = 8.72. 



Appendix A. The radius of gyration of the blobs: universal predictions 

In this appendix we wish to discuss the behavior of the radius of gyration of the blobs. If 
r g j($) is the radius of gyration of the i-th blob along the chain, the ratio r ffj j($)/_R 9 ($) 
is universal, being an adimensional ratio of large-scale properties of the polymer. It 
only depends on the position i of the blob along the chain, on the number n of blobs, 
and on the density through the polymer volume fraction $. Of course, this holds when 
the number of monomers L is large, otherwise scaling corrections should be taken into 
account. In general we have 



(A.l) 



where fi(n, and gi(n, are universal functions, A = 0.528(12), see Clisby, [5] is a 
universal exponent, and k a nonuniversal constant that does not depend on i, n, and $, 
but only on the model. In the polymer model we use at finite density, the Domb- Joyce 
model with w = 0.505838, the constant k is approximately zero, so that corrections 
decay with the next-to- leading exponent A2 ~ 1. 

An approximate expression for the n-dependence of the function fi(n, $) which 
works well for $ 1 is obtained as follows. In the large-L limit we have standard Flory 
scaling, R g = bL u and r gi = b'(L/n) u , with[5] v = 0.587597(7). Now assume that the 
blob shape and size is not influenced by the neighboring blobs, so that the size of the 
blob is equal to that of a free polymer with the same number of monomers. We can 
thus approximate b' pa b, so that r gi /R g = nT v . This formula is of course not exact, 
since blob-blob interactions cannot be neglected. Still, as we now show, it is reasonably 
accurate for $ 1. 

In order to compute r g ^/R g in the asymptotic limit, we determine Qi(L,n) = 
r g,i{ n )/Rg f° r L = Li = 600 and L = L 2 = 2400 in the Domb- Joyce model with 
w = 0.505838. Assuming corrections with exponent A 2 = 1.0(1), we estimate the 
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asymptotic (L — > oo) value as 

n ( \ L^Q l (L u n)-L^Q t (L 2 ,n) 

^i,as[n) — l/V.Z) 

L x - L 2 

The combination C(i, n, = n u Qi )aiS (n) for $ = is reported in figure [ATI as a function 
of (i — 1/2) /n for several values of n. Note that this quantity is always larger than 1, 
indicating that a blob of L/n monomers is more swollen than an isolated chain of the 
same degree of polymerization. This is due to the neighboring blobs which are entangled 
with the blob one is considering, causing swelling. Second, this effect is smaller for the 
boundary blobs since they have only one neighbor. The scaling r gi /Rg ~ n~ u holds 
quite well at zero density even for n = 4, with a proportionality constant which is only 
slightly larger than 1. In particular, for the boundary blobs we have f g ^jR g ~ LOln - ", 
while for the internal blobs f g ^/R g ~ 1.03n _I/ . If we average over all blobs and neglect 
end effects, we obtain relation (Q, which we used extensively in the text. 

The swelling effect is expected to increase as $ increases, since the higher the 
density the higher the blob-blob entanglement is. In figure IA~T1 we also report C(i, n, <3>) 
for $ = 8.72. There are here two notable differences with respect to the case $ = 0. First 
of all, end effects are small, indicating that much of the swelling is due to neighboring 
chains, consistently with the idea that for <3> > 1 polymers are strongly intertwined. 
Second, the n dependence of the scaling function fi(n, $) defined in Eq. ( 1A.ll) is not 
captured by the simple scaling form n~ v for our small values of n (of course, fi(n, $) 
scales as n~ v for n — > oo). 

Given the blob radii of gyration, using Eq. (@|, we can compute the ratio Rg^/Rg- 
For n = 4 we obtain 

' 0.89210(10) $ = 

R g>b (<S>) _ 0.88701(10) $ = 1.09 

R g (9) ] 0.87937(11) $ = 4.36 

k 0.8753(4) $ = 8.72 

Note that the $ dependence is tiny. At $ = 0, a good approximation for all n > 4 is 
given by 

Rg,b 



R g 



VI - kn- 2 » k = (1.03 - 0.04/n) , (A.4) 



which predicts for R g ^/ R g pa 0.8922 for n = 4, in good agreement with the result ( 1A.3I) . 



It is also interesting to consider the ratio r g (§,n)/r g (n), where r g is the average 
blob size in the asymptotic limit L — > oo (we perform the same extrapolation as done 
before for the ratios Qi). Results for several values of n and $ are shown in Table IA1 
and plotted in figure lA2l versus the blob volume fraction rjb = Cb/cl = 4nr g Cb/3. At least 
for rjb < 1 the data appear to depend only on r\b and to converge to 1 linearly as T]b —> 0: 
r g (§,n)/r g (n) pa 1 — 0.048^. Since r]b — > for n — > oo at fixed this result allows us 
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Table Al. Ratio r g (§,n)/r g (n) as a function of $ and n. 
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n= 10 


n = 20 


n = 30 


1.09 


0.982 


0.990 


0.994 


0.996 
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0.938 


0.962 


0.976 


0.982 
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Figure A2. Universal ratio r g ($,n)/r g (n) as a function of the blob volume fraction 
rjb- The dashed line is a linear fit of the data with rjb < 1. 



to predict the ratio the average of the Qi(n) defined above, as n — » oo. Indeed, 

we have 

r,(n) 4 # s (<&) ~ • J R 5 ($)' 1 • J 

The ratio R g (Q)/R g has been computed in several works. [Ml EE] For large $, R g /R g (<&) 
scales [IE] as 0.90$ a115 so that r g ($,n)/Rg(&) « 0.93n~^$ a115 . Note that scaling f[Q|) 
sets in for quite large values of n if $ is large. For instance, for $ = 8.72 it predicts 
n u Q = n u r g ($,n)/R g (&) w 1.23 for ri — >• oo, since [M] R g ($)/R g w 0.84 for this value 
of $. Hence, even for n = 30, see figure |A1~| we are still far from the asymptotic limit. 
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